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2 ! ABSTRACT 

I— I ' Context. The chemical analysis of secondary stars of low mass X-ray binaries provides an opportunity to study the formation processes 

(~| , of compact objects, either black holes or neutron stars. 

^ |. Aims. Following the discovery of overabundances of a-elements in the Keck I/HIRES spectrum of the secondary star of Nova Scorpii 

I ' 1994 (Israelian et al. 1999), we obtained VLT/UVES high-resolution spectroscopy with the aim of performing a detailed abundance 

O ' analysis of this secondary star 

, Methods. Using a ;^'2-minimization procedure and a grid of synthetic spectra, we derive the stellar parameters and atmospheric 

^2 ' abundances of O, Mg, Al, Ca, Ti, Fe and Ni, using a new UVES spectrum and the HIRES spectrum. 

, ^ , ' Results. The abundances of Al, Ca, Ti, Fe and Ni seem to be consistent with solar values, whereas Na, and especially O, Mg, Si and S 
are significantly enhanced in comparison with Galactic trends of these elements. A comparison with spherically and non-spherically 

^.f^ , symmetric supernova explosion models may provide stringent constraints to the model parameters as mass-cut and the explosion 

^ . energy, in particular from the relative abundances of Si, S, Ca, Ti, Fe and Ni. 

^.f^ • Conclusions. Most probably the black hole in this system formed in a hypernova explosion of a 30-35 M© progenitor star with a 

' mass-cut in the range 2-3.5 Mq . However, these models produce abundances of Al and Na almost ten times higher than the observed 

\^ \ values. 

Key words, black holes physics-stars:abundances-stars:individual (Nova Scorpii 1994, GRO J1655-40)-supernovae: general-X- 

, rays: binaries 
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1. Introduction cussions on the large amount of S, O, Mg and Si have consid- 

ered this object as a relic of a hypernova and Gamma Ray Burst 
The low mass X-ray binary Nova Scorpii 1994 (GRO J1655-40) ^^^^^^ 200O). More recendy, these abundance anomalies 



■ 

X 



IS one of the black hole transients that has attracted more atten- ^^^g compared with a variety of SN models, including stan 

H tion in the last decade. It was discovered on 1994 July 27 with ^^^^ ^^jj hypernova models (for several progenitor masses 

^_ BATSE on board the Compton Gamma Ray Observatory (Zhang geometries), and a simple model of the evolution of the bi- 

et al. 1994). The strong evidence that the compact object is a „^j.y ^^e pollution of the secondary (Podsiadlowski et al. 

black hole was presented by Bailyn et al. (1995) who suggested 2002). The observed anomalies were explained for a large range 

a mass function of f{M) - 3.16 + 0.15 Mq. A self-consistent j^^^^j parameters, providing substantial fallback and mixing 

analysis of the ellipsoidal light curves of the system has yielded between the fallback matter and the ejecta. 
masses of 5.4 + 0.3 and 1.45 + 0.35 Mq for the black hole and 
the companion star (Beer & Podsiadlowski 2002). 

Israelian et al. (1999) found evidence of the formation of the 

black hole in an explosive supernova (SN) event. Using a high- Mirabel et al. (2002) found additional independent evidence 

resolution spectrum obtained with the Keck I Telescope, they *at the black hole in this system formed in a SN explosion, 

found that the atmosphere of the secondary was enriched by a They determined the proper motion of the system using the 

factor of 6-10 in several a-process elements (O, Mg, Si and S). Telescope (HST) and computed its galactocentnc 

These elements, almost exclusively synthesized during super- J^^ ^y^tem moves with a runaway velocity of ~ 1 12 ± 18 

nova explosions, cannot be produced in a low-mass secondary s ' in a highly eccentric galactic orbit, 
star which must have been exposed to the supernova material 
ejected when the compact remnant was formed. Subsequent dis- 

Here we improve the chemical analysis of the secondary star 



Send ojfprint requests to: J. I. Gonzalez Hernandez including more elements like Na, Al, Ca, Ni and Li and refi- 

* Based on observations obtained with UVES at VLT Kueyen 8.2 m ning the abundance measurements of elements already studied 
telescope in programme 073.D-0473(A) by Israelian et al. (1999). 
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2. Observations 

We obtained 14 spectra of Nova Sco 94 with the UV- Visual 
Echelle Spectrograph (UVES) at the European Southern 
Observatory (ESO), Observatorio Cerro Paranal, using the 8.2 
m Very Large Telescope (VLT) on 2000 April 16 and 2000 June 
18, covering the specti-al regions ^^3300^500 A, ^^4800-5800 
A and ^^5800-6800 A at resolving power AI5A ~ 43, 000. The 
total exposure time was 5.6 hours. The spectra were reduced in a 
standard manner using the UVES reduction package within the 
MIDAS environment. In this paper, we will also make use of 
one spectrum taken with the HIRES spectrograph at Keck tele- 
scope, covering the spectral region /1/16800-8800 A, previously 
used by Israelian et al. (1999). This near-IR spectrum has been 
boxcar-smoothed in wavelength with steps of 0.3 A, providing a 
S/N -45-50 in the continuum. 



loop 
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Fig. 1. Radial velocities of Nova Scorpii 1994 folded on the orbital so- 
lution of the data with best fitting sinusoid. Individual velocity errors 
are < 1.5 km s"' and are not plotted because they are always smaller 
than the symbol size. The bottom panel shows the residuals of the fit. 



2.1. Revised orbital parameters 

We extracted the radial velocities by cross-correlating each 
UVES spectrum of the target with the spectrum of a F6IV 
template star (Shahbaz et al. 1999) taken from S'*N database 
(HIP 46853, Allende Prieto et al. 2004), using the software 
MOLLY developed by T. R. Marsh. A sine wave fit to the 
obtained velocities yields to the following orbital solution (see 
Fig.[B: r = -167.1 + 0.6 kms"', K2 = 226.1 ± 0.8 kms"', 
P = 2.62120 + 0.00014 d, and To = 2453110.5637 ± 0.0019 
d, where To is defined as the corresponding time of the closest 
inferior conjunction of the companion star, and the quoted un- 
certainties are at Icr. This orbital period, P, together with the 
velocity amplitude of the orbital motion of the secondary star, 
K2, leads to a mass function of f(M) - 3.16 + 0.03 Mq , con- 
sistent with previous results (Bailyn et al. 1995; Orosz & Bailyn 
1997). 

The derived radial velocity of the center of mass of the sys- 
tem disagrees somewhat (at the 3cr level) with previous studies 
(y = -150 + 19 km s-', Bailyn et al. 1995; y = -142.4 + 1.6 
km s ^ Orosz & Bailyn 1997; y = -I4I.9± 1.3 km s ', Shahbaz 
et al. 1999). Note that our high-resolution data has a factor 



greater than 5 better spectral resolution than the data these au- 
thors used. However, Shahbaz et al. (1999) estimate the uncer- 
tainty in their individual radial velocity measurements at typi- 
cally 6 km s"' , whereas our individual measurements have er- 
ror bars of roughly 1km s '. In addition, we have checked the 
accuracy of the wavelength calibration using the /16300.2 [Oi] 
sky line but unfortunately, its line profile appears distorted in 
each individual spectra of the target probably due to stellar lines, 
providing an estimate of its line center within 3 kms These 
uncertainties seem not to be enough to explain the difference be- 
tween both radial velocity measurements. 

2.2. Secondary spectrum 

The individual UVES spectra were corrected for radial veloc- 
ity and combined in order to improve the signal-to-noise ratio. 
After binning in wavelength with steps of 0.3 A the final UVES 
spectrum had a signal-to-noise ratio of 120-150 in the contin- 
uum in the red part of the spectrum. This spectrum is displayed 
in Fig. |2] Following Marsh et al. (1994), we computed the op- 
timal V sin / by subtracting broadened versions of a template 
star (HIP 46583), in steps of I km s and minimizing the resid- 
ual. We used an spherical rotational profile with linearized limb- 
darkening £ - 0.55 (Al-Naimiy 1978). The best fit corresponds 
to a vsini - 87 kms"', where the errors have been derived 
by assuming extreme cases for e = - 1 . These values are con- 
sistent with previous studies (vsin/ = 86^3 ^ kms~\ Shahbaz 
et al. 1999). Our derived rotational velocity, combined with our 
value of the velocity amplitude, K2, implies a binary mass ra- 
tio q = 0.329 + 0.047. Shahbaz (2003) might have determined 
a more accurate value of the mass ratio, q - 0.419 ± 0.028, by 
comparing synthetic spectra, which incorporate the secondary's 
Roche geometry, with an intermediate resolution spectrum of 
Nova Sco 94. 

2.3. Diffuse interstellar bands 

We obtained the spectrum of the interstellar medium (ISM) from 
our own observations. Due to the radial velocity of each indi- 
vidual exposure, in the average spectrum of all the individual 
spectra, the interstellar medium features are smoothed out. We 
subtracted the stellar features in each individual spectrum and 
combined them to obtain a spectrum that contains the ISM fea- 
tures. We tentatively identified in this spectrum 31 diff'use inter- 
stellar bands with equivalent widths typically in the range 10-65 
mA, although there are strong ISM features at 6060, 6202, 6283 
and 6613 A with equivalent widths of 148, 207, 555 and 175 
mA, respectively (Galazutdinov et al. 2000). In addition, the Na 
doublet at 5889 and 5895 A show equivalent width of 1.124 and 
1.006 A. We took care that these bands did not affect the lines 
selected for the chemical analysis described in the next section. 

3. Chemical Analysis 

3.1. Stellar Parameters 

In the normalized spectra of X-ray transients, although observed 
at quiescence, the stellar lines of the secondary star may be 
veiled by the flux coming from the accretion disc. This eff'ect is 
more relevant in secondary stars of low luminosity. In the opti- 
cal range, this veiling appears to drop toward longer wavelengths 
(e.g. Marsh et al. 1994). In our analysis, we tried to infer the stel- 
lar atmospheric parameters of the secondary star, using a tech- 
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Fig. 2. Observed spectrum of the secondary star of Nova Scorpii 1994 (top) and of a properly broadened template (HIP 46583, bottom). 



nique which compares a grid of synthetic spectra with the ob- 
served VLT/UVES spectrum, via a ;t'2-minimization procedure. 
This procedure takes into account the effect that veiling causes 
on the stellar lines and it has been applied in other X-ray bina- 
ries (A0620-00, Gonzalez Hernandez et al. 2004; Centaurus X- 
4, 2005; and XTE Jll 18-H480, 2006). For Nova Sco 94, we have 
also considered any possible veiling from the accretion disc. 

For this purpose, we firstly identified moderately strong and 
relatively unblended lines of the elements of interest in the high 
resolution solar atlas of Kurucz et al. (1984). We compute syn- 
thetic spectra with the local thermodynamic equilibrium (LTE) 
code MOOG (Sneden 1973), adopting the atomic line data from 
the Vienna Atomic Line Database (VALD, Piskunov et al. 1995) 
and using a grid of LTE model atmospheres (Kurucz 1993). 
We adjusted the oscillator strengths of the selected lines until 
reproducing both the solar atlas of Kurucz et al. (1984) with 
solar abundances (Grevesse et al. 1996) and the spectrum of 
Procyon with its derived abundances (Andrievsky et al. 1995). 
The changes we applied to the log gf values taken from VALD 
database are Alog gf < 0.2 dex. However, in some cases, we 
modified the log gf values even in the range 0.2-0.4 dex for 
some lines as for instance, the O i /16 156-8 A lines and the Mg i 
^5528 A and ^6318-9 A Hnes (see Tables |2] and O. 

We selected six spectral features containing in total 42 lines 
of Fe I with excitation potentials between 1 and 5 eV. We gen- 
erated a grid of synthetic spectra for these features varying as 
free parameters, the star effective temperature (Teff), the surface 
gravity (logg), the metallicity ([Fe/H]), and the veiling from the 
accretion disc, which was assumed to be a linear function of 
wavelength, and thus charaterized by two parameters, its value 
at 4500 A, /4500 = ^^df^c^/^se""' and the slope, otq. Note that 
the total flux is defined as Ftotai = f^disc + ^^.sec, where Fdhc and 



Table 1. Ranges and steps of model parameters 



Parameter 


Range 


Step 


Tss 


5500 -» 7500 K 


100 K 


log(g/cm S-) 


3. ^ 5 


0.1 


[Fe/H] 


-0.5 ^ 1 


0.05 


/4500 


0^ 0.6 


0.05 


mo 


^ -0.00027 


-0.00003 



Fsec are the flux contributions of the disc and the continuum 
of the secondary star, respectively. The microturbulence, ^, was 
fixed in each atmospheric model according to the calibration as a 
function of effective temperature and surface gravity reported by 
Allende Prieto et al. (2004). Such calibration has been derived 
for stars of the solar neighbourhood with similar metallicity as 
the secondary star in Nova Sco 94. For the stellar parameters of 
the secondary star, determined below, this calibration provides a 
microturbulence of ^ = 1 .9 km s" ' . 

The five free parameters were varying using the ranges and 
steps given in Table [T] In Fig. [3] we show the resulting his- 
tograms of these parameters corresponding to 1000 realizations 
of observed spectrum. We have not found any clear evidence 
of veiling from the accretion disc in the spectral region anal- 
ysed (/i/i5400-6800 A), as expected for the spectral type of 
the secondary star The comparison of this grid, using a boot- 
strap Monte-Carlo method, gives the most likely values Tetj = 
6100 ± 200 K, log(^/cm s^) = 3.7 + 0.2, [Fe/H] = -0.1 + 0.1, 
/4500 = 0.15 ± 0.05, and mo = -0.00012 + 0.00003. This 
metallicity determination is consistent with the values reported 
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Fig. 3. Distributions obtained for each parameter using Monte Carlo simulations. The labels at the top of each bin indicate the number of simula- 
tions consistent with the bin value. The total number of simulations was 1000. 



in the literature ([Fe/H] = 0.1 ± 0.2, Israelian et al. 1999; 
[Fe/H] = -0.1 ± 0.15, Buxton & Vennes 2001), whereas the 
effective temperature is slightly lower than the values obtained 
by these authors (T^ff = 6400 ± 250 K, IsraeUan et al. 1999; 
Teff = 6500 ±50 K, Buxton & Vennes 2001), although consistent 
within the error bars. We should remark that these authors used 
a different approach to derive the stellar parameters. Israelian 
et al. (1999) adopted the stellar parameters from previous cali- 
brations of the spectral type and luminosity classification of the 
secondary star as a F6III/F7IV star. On the other hand, Buxton 
& Vennes (2001) determined the effective temperature by fitting 
a lower quality spectrum, in particular, the wings of profile. 



3.2. Stellar abundances 

We inspected several spectral regions in both VLT/UVES and 
Keck I/HIRES spectra, searching for suitable lines for a de- 
tailed chemical analysis. Using the derived stellar parameters, 
we firstly determined the Fe abundance from each individual 
feature in the observed UVES spectrum (see Table|2]i. Since we 
knew of a potential enhancement of or-elements in the atmo- 
sphere of this star (Israelian et al. 1999), we computed new LTE 
model atmospheres with a-elements enhanced by H-0.4 dex with 
respect to solar using the Linux version (Sbordone et al. 2004) 



of the ATLAS code (Kurucz 1993). We used the new Opacity 
Distribution Functions (ODFs) of Castelli & Kurucz (2003) with 
the corresponding metallicity of the secondary star. In Fig. |4] 
we show some of the spectral regions analysed to obtain the Fe 
abundance. This figure also show the best synthetic spectral fit 
to the observed spectrum of a template star (HIP 46853 with 
reff= 6310 K, log §=3.91 and [Fe/H]= -0.08 dex) using the 
stellar parameters and abundances determined by Allende-Prieto 
et al. (2004). Thus, we only use as abundance indicators those 
features which were well reproduced in the template star. The 
chemical analysis is summarized in Table H] For the elements 
O, Mg, S, Si, Ti, previously studied by Israelian et al. (1999), 
we selected the most suitable lines in the UVES and HIRES 
spectra, whereas for other elements Na, Al, Ca, Ni, Fe and Li, 
we only used spectral features in the VLT/UVES spectrum. In 
Table |4] we provide the average abundance of each element ex- 
tracted from the analysis of several features. The number of fea- 
tures analysed for each element is also stated. Individual element 
abundances derived from each feature are given in Table |2] and 
Table [3] The errors on the element abundances show their sen- 
sitivity to the uncertainties on the effective temperature, Aj-^^, 
surface gravity, Aiog^, and the dispersion of the measurements 
from different spectral features, A^^. The errors A^- were esti- 
mated as Ao- = cr/ V^, where cr is the standard deviation of 
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Table 2. Adopted log gf values of important lines in the VLT/UVES spectrum. This line list do not included the weaker lines which 
these lines are blended with. The best fit abundance computed in LTE, [X/H], is also given for each feature. 



Species A{k) x i^^) ^og gf [X/H] Species ^ (A) ;r (eV) logg/ [X/H] 



Fei 


5462.970 


4.473 


-0.15 


- 


Si 


6748.682 


7.868 


-0.80 


- 


Fei 


5463.276 


4.434 


0.11 


-0.20 


Si 


6748.837 


7.868 


-0.60 


0.70 


Fei 


5466.396 


4.371 


-0.63 


- 


Si 


6757.007 


7.870 


-0.90 




Fen 


5534.847 


3.245 


-2.76 




Si 


6757.171 


7.870 


-0.31 


0.55 


Fei 


5535.418 


4.186 


-0.96 


0.15 


Si I 


5675.417 


5.619 


-1.23 


0.25 


Fei 


5554.895 


4.548 


-0.39 


-0.20 


Si I 


5675.756 


5.619 


-2.08 


- 


Fei 


5615.644 


3.332 


0.05 


-0.20 


Si I 


5684.484 


4.954 


-1.70 


0.60 


Fei 


6191.558 


2.433 


-1.52 


-0.10 


Si I 


6112.928 


5.616 


-2.25 


0.75 


Fei 


6230.723 


2.559 


-1.28 


-0.10 


Si I 


6113.120 


5.616 


-2.58 




Fei 


6393.610 


2.433 


-1.51 


-0.10 


Si I 


6131.565 


5.616 


-1.70 


0.25 


Fei 


6400.001 


3.602 


-0.41 


-0.20 


Si I 


6131.855 


5.616 


-1.70 


- 


Fei 


6592.914 


2.727 


-1.67 


-0.05 


Si I 


6142.204 


5.619 


-2.50 




Cai 


6102.723 


1.879 


-0.67 


-0.25 


Si I 


6142.483 


5.619 


-1.55 


0.50 


Cai 


6122.220 


1.886 


-0.01 


0.10 


Si I 


6145.016 


5.616 


-1.42 


- 


Cai 


6161.297 


2.523 


-1.35 


-0.25 


Si I 


6155.150 


5.619 


-0.82 


0.35 


Cai 


6162.180 


1.899 


0.37 


- 


Si I 


6155.693 


5.619 


-2.22 


- 


Cai 


6169.042 


2.523 


-0.78 


0.05 


Si I 


6237.319 


5.614 


-1.15 


0.70 


Cai 


6169.563 


2.526 


-0.45 


- 


Si I 


6238.286 


5.082 


-2.38 


- 


Cai 


6439.075 


2.526 


0.39 


0.00 


Si I 


6331.956 


5.082 


-2.40 


0.85 


Cai 


6471.662 


2.526 


-0.71 


0.15 


Si I 


6433.457 


5.964 


-1.62 


0.90 


Cai 


6475.236 


4.131 


-1.22 


-0.35 


Si I 


6721.848 


5.863 


-1.17 


0.70 


Cai 


6717.685 


2.709 


-0.62 


- 


Til 


5953.160 


1.887 


-0.52 


0.20 


Mgi 


5167.321 


2.709 


-1.03 


0.30 


Nal 


5682.633 


2.102 


-0.69 


0.70 


Mgi 


5172.684 


2.712 


-0.40 


0.00 


Nal 


5688.207 


2.104 


-0.25 


0.55 


Mgi 


5183.604 


2.717 


-0.18 


0.00 


Nal 


6154.226 


2.102 


-1.68 


0.05 


Mgi 


5528.405 


4.346 


-0.62 


0.60 


Nal 


6159.482 


2.104 


-1.36 


-0.05 


Mgi 


6318.717 


5.108 


-1.99 


0.30 


All 


6696.023 


3.143 


-1.83 


0.05 


iVlg 1 


6319 737 


5.108 


-2.25 




Al I 


6696 185 


4.022 


-1.83 




Mgi 


6319.495 


5.108 


-2.43 




All 


6696.788 


4.022 


-1.85 




Oi 


6155.966 


10.740 


-1.01 




All 


6698.673 


3.143 


-1.87 




Oi 


6156.756 


10.741 


-0.89 




Nil 


5475.425 


3.833 


-1.95 




Oi 


6156.776 


10.741 


-0.69 


0.90 


Nil 


5476.900 


1.826 


-0.89 


0.35 


Oi 


6158.176 


10.741 


-1.00 




Nil 


5625.312 


4.089 


-0.70 


0.10 


Oi 


6158.186 


10.741 


-0.41 




Nil 


5628.335 


4.089 


-1.27 


0.00 


Oi 


6300.230 


0.000 


-9.81 




Nil 


6116.199 


4.089 


-0.68 


0.45 


Si 


6743.440 


7.866 


-1.27 




Nil 


6116.199 


4.266 


-0.82 




Si 


6743.531 


7.866 


-0.92 


0.50 


Nil 


6118.094 


4.088 


-2.32 




Si 


6743.640 


7.866 


-1.03 




Nil 


6191.171 


1.676 


-2.55 


-0.10 


Si 


6748.573 


7.868 


-1.39 















Table 3. Same as Table|2]but for the important lines in the Keck I/HIRES spectrum. 



Species 


A (A) 


X(eV) 




[X/H] 


Species 


^(A) 


;r(eV) 


logg/ 


[X/H] 


Mgi 


7875.434 


5.932 


-2.13 




Oi 


7775.400 


9.110 


0.10 




Mg II 


7877.054 


9.996 


0.39 


0.75 


Oi 


8446.249 


9.521 


-0.46 




Mgi 


7877.478 


5.932 


-1.95 




Oi 


8446.359 


9.521 


0.23 


1.80 


Mgi 


7881.669 


5.933 


-1.51 


0.70 


Oi 


8446.759 


9.521 


0.01 




Mgi 


8712.676 


5.932 


-1.67 




Si 


8693.137 


7.870 


-1.37 




Mgi 


8712.689 


5.932 


-1.37 


0.60 


Si 


8693.931 


7.870 


-0.51 


0.90 


Mgi 


8717.803 


5.933 


-2.85 




Si 


8694.626 


7.870 


0.08 




Mgi 


8717.816 


5.933 


-1.85 




Si I 


8646.371 


6.206 


-1.94 




Mgi 


8717.825 


5.933 


-0.93 


0.65 


Si I 


8648.465 


6.206 


-0.30 


0.55 


Mgi 


8736.006 


5.946 


-1.93 




Si I 


8728.010 


6.181 


-0.61 


0.60 


Mgi 


8736.019 


5.946 


-0.69 


0.65 


Si I 


8728.594 


6.181 


-1.72 




Mgi 


8736.019 


5.946 


-1.97 




Si I 


8729.282 


6.181 


-2.30 




Mgi 


8736.029 


5.946 


-1.02 




Til 


8434.957 


0.848 


-0.88 


0.35 


Oi 


7771.960 


9.110 


0.45 


2.00 


Til 


8435.652 


0.836 


-1.02 




Oi 


7774.180 


9.110 


0.31 




Til 
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Fig. 4. Best synthetic spectral fits to the UVES spectrum of the secondary star in the Nova Scorpii 1994 system (bottom panel) and the same for a 
template star (properly broadened) shown for comparison (top panel). Synthetic spectra are computed for solar abundance ratios ([E/Fe] = with 
[Fe/H] = -0.1, dashed line) and best fit abundance (solid line). 



the measurements. The errors Aj-^i^ and Aiog^ were determined 

as Ar,, = (SlIi Ar,„,)/A^ and AiJgj = At,„,,j)/N. The 

Al and Li abundance were derived from one spectral line (see 
Fig. [8]l and the error associated to the dispersion of the mea- 
surements, cr, was assumed to be the average dispersion of Fe, 
Ca and Ni abundances, and in this case, Ao- - cr. The total er- 
ror given in Table |4] was derived using the following expresion; 

A,ot = ^/A^ + A?, -I- Ar . 

tot -y o- r^if logg 

3.2.1. Oxygen 

The atomic data for O i lines were adopted from Ecuvillon et al. 
(2006) where oscillator strengths were slightly modified in or- 
der to obtain a solar oxygen abundance of log ^(O)© - 8.74. In 
Fig. |5] we display one of the spectral regions analysed, where 
the forbidden line [Oi] /16300.2 A is present. Unfortunately, we 
found this line inconvenient to provide a reliable abundance: it 
is very weak at the effective temperature of the secondary star; it 
is blended with stronger lines of Si i /16299.6 A and Fez /J6301.5 
due to high rotational broadening of the star; and the strong 
narrow diffuse band (see the feature labelled DIB in Fig. |5]l at 
6283.8 A makes it difficult to place properly the continuum. 
An oxygen abundance of [O/H] ~ 0.9 (solid line in Fig. [Sj 
seems not to fit at all the [O i] feature which may be better repro- 
duced using a solar abundance synthetic spectrum (dashed line 
inFig.|5]). 

Thus, we determined the oxygen abundance from three fea- 
tures containing: Oi /16156-8 A lines, which appear severely 
blended with Si i, Ca i and Fe i lines, and the relatively isolated 



Oi /17771-5 A and Oi 8446 A triplets. In the left-bottom panel 
of Fig. |6]we show the Oi /16156-8 A lines in the VLT spec- 
trum of the target in comparison with several synthetic spec- 
tra in LTE. The best fit in LTE gives an oxygen abundance of 
[0/H]lte ~ 0.9. While LTE appears to be practically valid for 
the Oi ^6156-8 A (Takeda et al. 2005), the Oi ^7771-5 A and 
O I /18446 A suffers appreciable non-LTE effects (Ecuvillon et 
al. 2006). The NLTE correctional were kindly provided by N. 
Shchukina (private communication; Shchukina 1987; Ecuvillon 
et al. 2006). For the stellar parameters and oxygen abundance of 
the secondary star, NLTE corrections were estimated at -0.988 
and -0.926 for the near-IR Oi i7771-5 A and Oi /18446 A 
triplets, respectively. 

3.2.2. Magnesium 

The abundance of magnesium was determined from 5 Mg i spec- 
tral features in the optical VLT/UVES spectrum, including the 
three features of the Mgib triplet at 5167-83 A, 4 Mgi spec- 
tral features in the near-IR Keck I/HIRES spectrum, where Mg i 
lines are isolated or weakly blended with other lines, and a Mg ii 
spectral feature at 7877 A. In Fig.|5]and the right-bottom panel of 
Fig.|6] we show some of the Mg i features included in the abun- 
dance analysis. Note that the best fit synthetic spectra displayed 
in these figures show the best fit abundance of each feature in- 
stead of the average abundance of each element given in Table|4] 
We should remark that the location of the continuum near the 
near-IR lines, where the S/N~ 50, is less accurate, than that of 



' AnlTE = log f(X)NLTE - log e(X)LTE 
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Fig. 5. The same as in Fig. 4, but for the spectral range /1/16275-6400 A. 
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Fig. 6. The same as in Fig. 4, but for the spectral range /i/i6141-6172 A (left panels) and /1/15521-5559 A (right panels). 



the optical lines, where the S/N~ 150. Thus, we estimate an av- 
erage uncertainty of ~ 0.27 dex in the Mg abundance due to the 
uncertainty on the continuum location in the near-IR spectrum, 
whereas an uncertainty of ~ 0.05 - 0.1 dex is expected from the 
optical features. This might explain the high dispersion of the 
Mg measurements given in Table |4] On the other hand, the near- 
IR Mg features are less sensitive (< 0.1 dex) to the errors on the 
stellar parametes than the optical lines, specially the Mg ib triplet 
whose errors related to the uncertainties on the stellar parameters 

are At^ ~ 0.2 - 0.3 dex, and Aiog s, 0.1 dex. More details on 

the possible uncertainties on the analysis of Mg ib lines are given 
in § 13.3.31 In addition, some of the Mg lines analysed are known 
to exhibit non-LTE effects. NLTE abundance corrections, Anlte, 
for Mgl lines have been studied by Zhao et al. (1998, 2000) in 



the Sun. Among the lines we have analysed, they only reported 
corrections for the Mgi /15528 A and the Mgib /15 167 - 5183 
A, which are expected to show NLTE abundance corrections of 
roughly — 1-0.05 dex in the Sun and, in general, less than 0.1 dex 
in dwarfs and subdwarfs (Shimanskaya et al. 2000). On the con- 
trary, Abia & Mashonkina (2004) estimated NLTE corrections 
of -0.02 < Anlte < -0.07 in the Sun for the Mg ii .17877 A. 
Although these NLTE corrections are small, they should be in- 
vestigated for the case of the abundance and stellar parameters 
of the secondary star in this system. 
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Table 4. Chemical abundances of the secondary star in 
Nova Scorpii 1994 



Element 


[X/H]i^TE 


(T 


K 






A,o, 


lines 


0"= 


0.91 


0.07 


0.04 





0.07 


0.09 


3 


Na 


0.31 


0.37 


0.18 


0.19 





0.26 


4 


Mg 


0.48 


0.29 


0.09 


0.10 


-0.03 


0.15 


10 


Al 


0.05 


0.15 


0.15 


0.1 





0.18 


1 


Si 


0.58 


0.21 


0.06 


0.05 


0.01 


0.08 


12 


s 


0.66 


0.18 


0.09 


-0.04 


0.06 


0.12 


4 


Ca 


-0.02 


0.17 


0.07 


0.12 


0.01 


0.14 


7 


Ti 


0.27 


0.11 


0.07 


0.2 


0.05 


0.22 


2 


Fe 


-0.11 


0.11 


0.03 


0.09 


0.01 


0.10 


9 


Ni 





0.10 


0.06 


0.2 


0.05 


0.21 


5 


Lif 


2.13 


0.15 


0.15 


0.15 





0.21 


1 



" The solar element abundances were adopted from Grevesse et al. 
(1996) and Ecuvillon et al. (2006): log 6(0)o = 8.74, log e(Na)o = 6.33, 
log6(Mg)e = 7.58, loge(Al)o = 6.47, log6(Si)o = 7.55, loge(S)o = 
7.21, loge(Ca)o = 6.36, log6(Ti)o = 4.99, log6(Fe)Q = 7.50, 
log6(Ni)o = 6.25 

Element abundances (calculated assuming LTE) are [X/H] = 
log[A'(X)/yV(H)]„„ - log[A'(X)/yV(H)]su„, where A'(X) is the number 
density of atoms. Uncertainties, Atot, are at Icr level and take into ac- 
count the error in the stellar parameters and the dispersion of the mea- 
surements. 

The errors from the dispersion of the best fits to different features, see 
text. 

Number of features analysed for each element. 
NLTE abundance of oxygen, see text. 
' Li abundance is expressed as: 

l0ge(Li)NLTE = log[A'(Li)/A'(H)]NLTE + 12 



3.2.3. Silicon 

The Si abundance was obtained from the 12 features with Si I 
lines often blended with other lines, some of them are displayed 
in Figs. |7] and [8] The dispersion of the abundance measurements 
from different lines is 0.2 dex, although the average abundance 
has a good accuracy due to the number of features analysed. As 
shown in Fig.|5] we found difficult to reproduce Sin lines, sur- 
prisingly strong in comparison with Si i lines, with the average 
Si abundance, so that, we decided not to consider them as abun- 
dance indicators. Shchukina & Trujillo Bueno (2001) pointed 



out that NLTE effects for Si (with an ionization potential of 8.15 
eV, slightly higher than that of Fe, 7.87 eV) would likely be sim- 
ilar to those found for Fe, which suffers from over-ionization, in 
a three-dimensional solar photospheric model. In addition, for 
species where over-ionization is prominent, the NLTE effects 
normally increase toward higher T^ff and lower log g (Asplund 
2005). 

Wedemeyer (2001) has performed a NLTE analysis of Si in 
the Sun and argued that most of the levels of Si i show almost 
negligible deviations from LTE, whereas Sin levels are over- 
populated with respect to LTE. In particular, this author pre- 
dicts NLTE corrections of -0.01 dex for one of the Sii lines we 
have analysed at 5684 A. For the Si ii hnes at 6347 A and 6371 
A the NLTE corrections are -0.097 dex and -0.064 dex for the 
Sun, respectively, and therefore 6 to 9 times higher for the Si ii 
lines than for the Sii lines. These NLTE effects might be sig- 
nificantly larger for the Si ii lines in the secondary star, giving 
rise to large NLTE corrections. The synthetic profile computed 
in NLTE would make stronger Si ii lines, thus providing a better 
description of the observations. However, this conclusion must 
be clarified in forthcoming investigations. Note that, in fact, the 
synthetic spectral fits reproduce quite well the Si ii features in the 
template, but these features are much stronger in the spectrum of 
the secondary star in Nova Sco 94. 



3.2.4. Sulpliur 

In Fig. [8] we show another spectral region with the strong S i 
/16743-57 A lines. The overabundance of S is clearly seen in 
the spectrum of Nova Sco 94 when compared with the template. 
These three strong and unblended lines of sulphur clearly con- 
firm the overabundance of sulphur in the secondary star of this 
system. We also included one S i feature at 8694 A, found in 
the Keck spectrum, in the S abundance determination. Takeda- 
Hidai et al. (2002) and Takeda et al. (2005) have reported neg- 
ligible corrections for S i /i8694 A lines. For the stellar param- 
eters, metallicity and abundance ratio [S/Fe]~ +0.6 dex of the 
secondary star, NLTE coiTections are comprised in the range 
-0.04 < Anlte < -0.08 dex. These coiTections were not con- 
sidered in the sulphur abundance given in Table H] since they are 
almost negligible compared to the dispersion of the LTE abun- 
dances derived from different S i lines. In addition, NLTE cor- 
rections for the S i /16743-57 A lines are not available in the 
literature. 
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Fig. 8. The same as in Fig. 4, but for the spectral range /l/?6690-6760 A. 



3.2.5. Calcium 

Calcium is one the a-elements that could not be measured in 
the analysis of the Keck spectrum reported in Israelian et al. 
(1999). The abundance of calcium was obtained from the anal- 
ysis of 7 almost unblended and strong spectral features in the 
VLT spectrum. Some of these features are displayed in the left- 
bottom panel of Figs. |6l|7]and[8] NLTE abundance corrections 
have been reported for the Ca i /16102, ^6169, .16439 and .16471 
A lines (Mashonkina et al. 2007), being -0.05, -0.01, -0.13, - 
0.06 dex, respectively for a model of reff= 6000 K log g= 4 
dex and [Fe/H]= dex. This results on an average correction 
Anlte < -0.05 dex which does not have any significant impact 
on the Ca abundance. 

3.2.6. Titanium 

Israelian et al. (1999) could not measure a reliable abundance of 
Ti using the Keck spectra since they did not find spectral features 
with relatively unblended Ti i lines. Unfortunately, the signal-to- 
noise ratio in the spectral region from 5920 to 5960 A of the 
VLT/UVES spectrum of Nova Sco 94 where there are some Ti 
features is not enough to provide reliable abundances. This spec- 
tral region is close the edge of one of the arms of the UVES 
spectrograph. Thus, we finally selected one of these features in 
the VLT spectrum and another relatively unblended feature of 
Ti in the Keck spectrum. We emphasize that the Ti abundance 
must be considered with caution, although the UVES spectrum 
seems to indicate that Ti abundance is not as enhanced as other 
a-elements like O, Mg, Si and S, but it shows an abundance al- 
most consistent with solar if we take into account the error bars 
(see TablelH. 



3.2.7. Tine odd-Z elements: Na, Al 

The All /16696-8 A lines appear to be surprisingly weak 
since one could expect to obtain an overabundance of alu- 
minum according to supernova yields (Umeda & Nomoto, 2002, 
2005; Tominaga, Umeda & Nomoto 2007; Maeda et al. 2002). 
Baumuller & Gehren (1996, 1997) show that this doublet is 
not sensitive to NLTE conditions, with expected corrections 
Anlte < +0.05 for models with similar stellar parameters as the 
secondary star. 

The sodium abundance was determined from 4 features, typ- 
ically blended with other lines (see Table |4|. Its abundance is 
not as high as expected from yields of SN models which sug- 
gest a similar behaviour as aluminum. Baumuller et al. (1998) 
predict NLTE corrections Anlte < -0.03 and Anlte < -0.05 
for the Al i /15682-8 and /16 154-60 lines, respectively, which are 
not relevant in the abundance analysis. 

3.2.8. Lithium 

The best fit to the Li i /16708 A feature provides an LTE abun- 
dance of log e(Li)LTE - 2.16 + 0.21. We estimated the non-LTE 
abundance correction for this element, from the theoretical LTE 
and non-LTE curves of growth in Pavlenko & Magazzil (1996). 
We found Anlte - -0.03. Due to the weakness of the absorp- 
tion we consider this abundance estimate given in Table |2] as an 
upper limiB 

According to the effective temperature of the secondary star 
(Tes - 6100 + 200 K), slightly far from the lithium gap (between 
6400 - 6800 K, Boesgaard & Tripicco 1986), this Li abundance 

^ The template star has log6(Li)LTE ~ 
3.3, whichdoesnotaffectourre suits 
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Fig. 9. Upper panel: Observed UVES spectrum of the secondary star in comparison with the average synthetic spectrum computed with LinBrod 
(dashed line) for all orbital phases and the synthetic spectrum computed with MOOG (solid line). An abundance of [Mg/H]~ 0.3 dex has been 
adopted for all synthetic spectra. Bottom panel: Synthetic spectra computed with LinBrod for different orbital phases. 



is consistent with that of main-sequence F type disc stars with 
ages in the range ~ 1 - 3 x 10^ yr with the similar metallicity 
(e.g. Boesgaard & Tripicco 1987; Balachandran 1991). 

3.3. Comparison with the results of Foellmi et al. 

In a recent paper, Foellmi, Dall & Depagne (2007) compared 
the results on the element abundances reported by Israelian et 
al. (1999) using the UVES spectra reported in this paper, taken 
in 2004, and additional VLT/UVES spectra taken in 2006. They 
confirm an enhanced oxygen abundance and a relatively low Ti 
and Ca abundance in the secondary star which is consistent with 
our results, but also claim that other a-elements do not show 
over-abundances in clear disagreement with our findings for Mg, 
Si and S. Their conclusions are based on equivalent widths of 
different spectral features as well as the comparison of synthetic 
spectra to the observed Mg ib triplet at 5167-83 A, the S i triplet 
at 9228 A and some other features of Mg i. Si i and S i in the 
spectral region from 8670 to 8740 A. 

We have revised the equivalent widths presented in Table 1 
of Foellmi et al. (2007) in our VLT/UVES spectra obtained in 
2004 and found slightly larger values which might be related to 
the continuum location. In addition, we find it confusing the sys- 
tem velocity of 88 km s"' reported by Foellmi et al. (2006) 

which they also used in Foellmi et al. (2007), in comparison 
with our value of ~ -167 kms ' which is comparable to val- 
ues previously reported in the literature. However, they obtained 
an orbital semiamplitude of K2 = 225 + 10 km s"' which is less 
accurate than our value, K2 = 226.1 ± 0.8 km s"'. Their larger 
error on the orbital semiamplitude could give rise to large errors 



on the radial velocities applied to the individual spectra which 
might have slightly smoothed out the stellar lines in the average 
spectrum of the secondary star. This might explain why these 
authors found a rotational velocity of v sin / ~ 94 km s"' slightly 
larger than our value of v sin / ~ 87 km s ' . 

3.3.1. Sulphur 

Foellmi et al. (2007) claim that sulphur is not enhanced in 
the secondary star Indeed, we demonstrate in our paper that 
the sulphur features at 6743-57 A are very well detected with 
S/N~ 150 (as it is shown in Fig.[8]of this paper) and lead to a 
large overabundance of this element. The new features we are 
referring in our new paper basically confirm the sulphur over- 
abundance in a completely independent way to Israelian et al. 
(1999). Foellmi et al. did not note these sulphur features in their 
paper which should be certainly present in the spectra they used. 
In addition, in Figs. 1 and 3 in Foellmi et al. (2007), the loca- 
tion of the continuum appears improperly low. The S i features 
at 8693 and 9228 A would lead to a larger abundance if the con- 
tinuum is suitably placed. We should remark that due to the large 
rotational broadening of the secondary star, roughly 87 km s"', 
slight changes of the continuum location could give rise to con- 
siderably differences in the derived element abundances. The in- 
spection of the lower panel of Fig. 1 in Foellmi et al. (2007) 
shows that the Fei features, which are not expected to be en- 
hanced, are not well reproduced by the synthetic spectrum be- 
cause the continuum has been placed too low. After correction 
of the continuum in this figure, the S i and Si i features at 8680 
and 8694 A would require higher abundances. 
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3.3.2. Silicon 

Foellmi et al. (2007) also argue that silicon is not enhanced. We 
should note that we used 12 features of Si i in our analysis, and 
10 of these features are located in regions of the UVES spec- 
trum were the signal-to-noise ratio is roughly 150. We should 
also remark that Mg, Si and S show an abundance dispersion 
from the measurement of different features of ~ 0.2 - 0.3 dex, 
which indicates that there are several features which provide 
lower abundances than the average which is given in Table |4] 
It is therefore convenient to derive the element abundances from 
different features, and then provide the mean and the dispersion. 
Foellmi et al. (2007) have just compared one synthetic spectrum 
with the observations without treating each feature individually. 
It is therefore not straightforward to work out the reasons for 
the differences. We should remark that we have carried out a 
differential analysis with respect to the template star and the 
sun. Thus our atomic data have been checked by comparing the 
solar atlas with a synthetic spectrum computed with the same 
code and adopting the solar atmospheric parameters and abun- 
dances. Foellmi et al. do not provide information on their atomic 
data. Furthermore, we treated carefully the normalization of the 
observed spectrum in comparison with the spectrum of a tem- 
plate star with similar stellar parameters and known abundances. 
Thus, we only used as abundance indicators those features which 
were well reproduced in the template star. Moreover, it is also 
clear there are some aspects that are difficult to explain in the 
data as the strong Si ii lines in the optical spectrum (as it is shown 
in Fig.|5]of this paper), which would require an overabundance 
of [Si/H]~ 1 .75 dex in LTE. We noted but did not use these fea- 
tures in our analysis which would add for overabundance. We 



should note that decreasing the surface gravity by 0.2 dex, which 
might be expected from the size of the Roche lobe (see § |3.3.3t . 
does not improve significantly the agreement between the syn- 
thetic and the observed Sin lines. In addition, the Fen /16147 
A feature is not either well reproduced using this value of the 
surface gravity, despite that Fe ii lines are not expected to show 
deviations from LTE. 



3.3.3. Magnesium 

Finally, Foellmi et al. (2007) also discussed the apparently weak- 
ness of the Mgib triplet at 5167-83 A. They infer a Mg abun- 
dance in LTE of [Mg/H]~ 0.1 dex. Our spectroscopic determi- 
nation of the abundance provided by the Mg ib triplet using the 
code MOOG is consistent with this result (see Table|2]i. As dis- 
cussed in § 13.2.21 the Mg ib is quite sensitive to possible errors 
on the stellar parameters, so we prefer not to rely on a magne- 
sium determination based just on this feature. We show in what 
follows that the Mg Ib triplet is also sensitive to deviations from 
sphericity in the stellar atmosphere. The secondary star is fill- 
ing its Roche lobe, and therefore is elongated towards the in- 
ner Lagrange point from which the star loses matter onto the 
black hole. Thus, the temperature, surface gravity and Doppler 
shift at line formation are not constant over the whole surface of 
the star, and consequently, the profile of certain lines depends 
on the orbital phase. In order to see if this is the case for a 
given spectral feature we have used the program LinBrod which 
computes synthetic spectra for each orbital phase taking into ac- 
count the Roche lobe symmetry of the stellar surface (Bitner & 
Robinson 2006). The input parameters of this program are the 
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Fig. 11. Abundance ratios of the secondary star in Nova Scorpii 1994 (blue wide cross) in comparison with the abundances of G and K metal-rich 
dwarf stars. Galactic trends were talsen from Ecuvillon et al. (2004), Ecuvillon et al. (2006) and Gilli et al. (2006). The size of the cross indicates 
the 1-sigma uncertainty. Filled and empty circles correspond to abundances for planet host stars and stars without known planet companions, 
respectively. For the abundance of oxygen in metal-rich dwarfs we have only considered abundance measurements in NLTE for the triplet O i 
7771-5 A. The dashed-dotted lines indicate solar abundance values. 



average effective temperature of the star, which was assumed to 
be our spectroscopic estimate of the effective temperature, and 
the orbital parameters. We also adopted an orbital inclination of 
~ 69° (Beer & Podsiadlowski 2002) and mass ratio of q ~ 0.419 
(Shahbaz 2003). The code LinBrod divides the stellar surface in 
tiles, for which a local temperature and gravity is derived, uses 
the code MOOG to compute a synthetic spectrum for each tile, 
and finally integrates over the whole stellar surface. The average 
profile of Mg ib lines over all the observed orbital phases com- 
puted with LinBrod is weaker than the spherical profile com- 
puted with the code MOOG (see Fig.|9j. The main reason of this 
difference is that the average surface gravity derived by LinBrod, 
which is shared by almost the 80 per cent of the tiles and only 
depends on the dynamical parameters, is logg - 3.5 + 0.1, being 
lower that our estimate (log g - 3.7 + 0.2) but marginally consis- 
tent within the error bars. The Mg ib lines at 5172 A and 5183 A 
become weaker for lower values of log g. The Mg ib line at 5 167 
A is blended with lines of other elements that have different sen- 
sitivity to changes of the stellar parameters and they somehow 
compensate the sensitivity of the Mgib 5167 A line. This figure 
also shows how the profile depends on the orbital phase, mostly 
due to the different velocity profile provided by the asymmetry 
of the stellar surface seen by the observer. 

In Fig. |9] we adopted an abundance of [Mg/H]~ 0.3 dex 
which is almost able to reproduce the observed features when 
the code LinBrod is used. Note that we have taken into account 



a veiling factor of 7% estimated at this wavelength according to 
our determination of stellar and veiling parameters (see § 13. It . 
This abundance would be in agreement with the Mg i feature at 
6319 A which provides [Mg/H]~ 0.3 using the code MOOG and 
is less sensitive to the surface gravity (see Table|2|. 



In Fig. [To] we displayed another spectral region where the S i 
6743-57 A are located. In this case, in contrast with the Mg ib 
lines, the synthetic spectrum computed with LinBrod is practi- 
cally similar to that computed with MOOG, which indicates that 
these S i lines are not affected by the asymmetry of the stellar sur- 
face because they are hardly sensitive to the surface gravity. The 
Mg features analysed in the near-IR Keck spectrum are equally 
not affected because they are not sensitive to small variations of 
the stellar parameters and in addition, this spectrum was taken at 
an orbital phase close to 0, when the star looks almost spherical. 



We should remark that in the average Mg abundance given in 
Table m we do not applied any correction to the Mg abundances 
derived from the Mg ib lines using the code MOOG, which are 
given in Table|2] The use of the program LinBrod would provide 
a slightly higher average Mg abundance and slightly smaller dis- 
persion. 
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Fig. 12. Left panel: Observed abundances (filled circles with the error bars) in comparison with the expected abundances in the secondary star after 
having captured the 70% of the matter ejected within the solid angle subtended by the secondary from a spherically symmetric supernova explosion 
of kinetic energy = 10^' ergs for two different mass-cuts, M^u, = 2 Mq (solid line with open circles), and Mem = 2.54 Mg (dashed-dotted line 
with open circles). Right panel: the same as left panel but for a spherically symmetric hypernova explosion of kinetic energy Ej^ =2Qx 10^' ergs 
for two different mass-cuts, M^m = 3.04 Mg (solid line with open circles), and M^„^ = 4.05 Mq (dashed-dotted line with open circles). 



4. Discussion 

4.1. Heavy elements 

The abundances ratios of O, Mg, S, Si and Na with respect to 
Fe in the secondary star of Nova Sco 94 are considerably higher 
than those in stars of similar Fe content (see Fig. fTTT i while, on 
the contrary, Al, Ca, Ti and Ni seem to be consistent with the 
Galactic trends. As can be seen in Table |4] the uncertainties 
induced by effective temperature and gravity are considerably 
diminished when dealing with abundance ratios and the major 
source of error in [X/Fe] is associated with the dispersion, Aq-, 
of abundances obtained from different features of the same ele- 
ment. Thus, the error bars of the abundance ratios displayed in 
Fig, nn were estimated as: 

(A[X/Fe])2 = A^ X + (A7-„„x - A7-„„Fe)' + 
(Alogg,X - Aiogg,Fe)^ + A^ pg 

Podsiadlowski et al. (2002) showed that a relatively simple 
model of pollution of the secondary star from a supernova ex- 
plosion can explain most of the observed (T-element abundance 
enhancements in the secondary star 

Since we have extended the number of elements analysed in 
the secondary star, we will intend to extract more information on 
the SN explosion by comparing our new results with yields from 
current supernova models kindly provided by Umeda, Nomoto 
et al. (2002, 2005). 

4.1.1. Spherical SN explosion 

We have considered spherically symmetric supernova and hy- 
pernova models for a 30 M© progenitor (8.5 He core, Umeda & 
Nomoto, 2002, 2005; Tominaga, Umeda & Nomoto 2007) and a 
~ 2 Mq secondary star, initially placed at an orbital distance of 
~ 10 Rq (after the binary orbit has been recircularized), which 
needs to capture almost 70% of the matter ejected within the 
solid angle subtended by the secondary as seen from the helium 
star to achieve the observed abundances. We adopted an initial 
black hole mass of Mbh,! ~ 5 Mq , before the secondary star 
begins to transfer matter onto the compact object. Note that the 



dynamical estimate of the black hole mass provides Mbh,! ~ 5.4 
Mq (Beer & Podsiadlowski 2002). We have also assumed differ- 
ent mass-cuts (M^ut, the mass that initially collapsed forming the 
compact remnant) and all the fallback material (Mfau, amount 
of mass which is eventually accreted by the compact core) to 
be completely mixed with the ejecta. We have performed sim- 
ulations of the expected abundances of the secondary star after 
capturing matter from the SN ejecta as in Gonzalez Hernandez et 
al. (2004). In our simulations we varied the mass-cut from 1.35 
to 5.05 Mq in steps of ~ 0.5 and 1 M© (depending on the mass 
bins given in the explosion models), and all the simulations keep 
Mbh.i - A^cut + A/faii. In addition, all the matter captured is com- 
pletely mixed with the whole material of the secondary star. In 
Fig. [12] we show the observed abundances in comparison with 
the expected abundances for two different mass-cuts. The error 
bars in Fig. [T2l correspond to the total errors given in Table |4] 
Note that some of the elements are highly sensitive to mass- 
cuts, for instance, the supernova model (left panel in Fig. [12] ) 
predicts a strong dependence of Ni, Ca and S abundances on 
the selected mass-cuts since the explosive burning is only able 
to create these elements at relatively deep layers in the explo- 
sion. It would be acceptable to find an agreement between the 
observed abundances and one of the two models displayed or at 
least to find the observed values in a middle location between 
both predictions. This comparison using both spherically sym- 
metric supernova and hypernova models allow us to extract the 
following conclusions: 

a) The relatively low abundances of Ca, Ti, Ni and Fe may in- 
dicate that the mass-cut was in the range 2-3.5 Mq for ex- 
plosion energies (1 - 20) x 10^' ergs. 

b) These explosion models suggest that the higher the energy 
the larger amounts of elements like S, Si, Ti and Ca, and even 
Fe and Ni are created. Thus, the substantial enhancement of 
Si and S in the secondary star favours higher explosion ener- 
gies and hence, the hypernova scenario (Brown et al. 2000). 

c) In addition, even for higher energies, the amount of S de- 
creases significantly above a mass-cut of 3.5 M© confirming 
the need of efficient mixing processes between the fallback 
matter and the ejecta (Kifonidis et al. 2000). 
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Fig. 13. Left panel: Observed abundances (filled circles with the error bars) in comparison with the expected abundances in the secondary 
star after having captured the 10% of the matter ejected within the solid angle subtended by the secondary from a non-spherically symmetric 
supernova explosion of kinetic energy Eic = 10^^ ergs for two different mass-cuts, M^„t = 2.41 Mq (solid line with open circles), and M^„t = 2.66 
Mq (dashed-dotted line with open circles). This model corresponds to the matter ejected in the equatorial plane of the primary where we assumed 
that the secondary star is located (see more details in Gonzalez Hernandez et al. 2005b). Right panel: the same as left panel but in this model we 
have assumed complete lateral mixing where all the material within given velocity bins is assumed to be completely mixed. Two simulations are 
shown for two different mass-cuts, M^ut = 2.41 M© (solid line with open circles), and M^h = 5.14 Mq (dashed-dotted line with open circles). 



d) The abundances of O, Mg, Al and Na hardly depend on the 
explosion energy and mass-cut, but interestingly Al and Na 
must be more overabundant than O and Mg in contrast with 
the observed abundances. 

4.1.2. Aspherical SN explosion 

The derived heliocentric radial velocity of ~ -167.1 kms ' is 
much larger than the system kick (~ 50 km s"') gained from a 
spherically symmetric supernova explosion. Brandt et al. (1995) 
suggested the black hole in Nova Sco 94 could have formed in a 
two-stage process where the collapse first leads to the formation 
of a neutron star accompanied by a supernova kick, which sub- 
sequently converted into a black hole by fallback. Podsiadloski 
et al. (2002) also investigated the expected abundances in the 
companion star after pollution of nucleosynthetic products in a 
non-spherically symmetric supernova. They found good fits to 
the observed abundances reported in Israelian et al. (1999). In 
their simulations they also took into account possible losses of 
mass of the He core progenitor in a stellar wind before the explo- 
sion and the pollution because of the capture of wind material. 
This allowed them to use He core progenitor with masses as high 
as 16 Mq , hence reducing the matter ejected in the SN explosion 
and preventing from the disruption of the system. 

We have also inspected an aspherical explosion model of a 
16 Mq He core progenitor with an explosion energy of 10 x 10^' 
ergs from Maeda et al. (2002). We did not consider mass losses 
from the helium star via stellar winds and hence the greater 
amount of ejected matter only requires to capture roughly 10% 
of the matter ejected within the solid angle subtended by the sec- 
ondary. We computed the expected abundances in the secondary 
star after pollution from the SN products assuming that the sec- 
ondary star was located in the equatorial plane of the helium star 
before the explosion. In this case, the mass-cut varies from 2.41 
to 5.14 Mq in steps of ~ 0.6 Mq . Fig. [13] shows the compari- 
son of the observed abundances with the expected abundances 
in the secondary star after pollution from the SN ejecta. The er- 
ror bars in Fig.[T3]coiTespond to the total eiTors given in Table|4] 



As in the spherical case we show two models and two different 
mass-cuts. The left panel reflects the composition of the material 
ejected in the equatorial plane while in the right panel we have 
considered complete lateral mixing (Podsialowski et al. 2002), 
i.e. the ejected matter is completely mixed within each velocity 
bin. Note that in the left panel we selected two close values of 
the mass-cut while in the right panel the whole dependence on 
the mass-cut is shown. The model with complete lateral mixing 
tends to enhance all the element abundances for each mass-cut. 
The analysis of the aspherical case provides the following con- 
clusions: 

a) The aspherical model predicts too low abundances of S for 
any mass-cut unless complete lateral mixing is assumed. 

b) However, in the model with complete lateral mixing, in spite 
of the fact that Na abundance decreases until being almost 
marginally consistent with the observed value, the observed 
Al still remains too low in comparison with its expected 
abundance. 

c) Moreover, this model increases Ca, Ti, Fe and Ni abundances 
by a factor of 2-3 above the observed values. 

A complete grid of models with a wide range of progenitor 
masses and explosion energies may be required to study in detail 
the abundance pattern of the secondary star in this system. 

5. Conclusions 

We have presented a VLT/UVES high resolution spectroscopy 
of the black hole binary Nova Sco 94. The individual spectra 
of the system allowed us to derive an orbital period of P = 
2.62120 ± 0.00014 days and a radial velocity semiamplitude of 
the secondary star of K2 - 226.1 + 0.8 km s The implied up- 
dated mass function is /(M) = 3.16 + 0.03 Mq , consistent with 
previous values reported in the literature. The inspection of the 
high-quality averaged spectrum of the secondary star provides a 
rotational velocity of v sin / = 87^^ km s"', and hence a a bi- 
nary mass ratio q = 0.329 ± 0.047. The derived radial velocity. 



J. I. Gonzalez Hernandez, R. Rebolo, and G. Israelian: The Black Hole Binary Nova Scorpii 1994 (GRO J1655-40): 



15 



y = -167.1 ± 0.6 km s" , of the center of mass of the system 
disagrees, at the 3<t level, with previous studies. 

We have performed a detailed chemical analysis of the sec- 
ondary star. We applied a technique that provides a determina- 
tion of the stellar parameters, taking into account any possible 
veiling from the accretion disc. We find Tefr = 6100 ± 200 K, 
log(^/cms2) = 3.7 ±0.2, [Fe/H] = -0.1 ±0.1, and a disc veiling 
(defined as Fdisc/^totai) of less than 10% at 5000 A and decreas- 
ing toward longer wavelengths. 

We have revised the chemical abundances of O, Mg, S, Si, Ti 
and Fe already reported in Israelian et al. (1999) and determined 
new element abundances of Na, Al, Ca, Ni and Li. The element 
abundances are typically higher than solar except for Fe, Ca and 
Ni, and in some cases significantly enhanced (e.g. O, Mg, S and 
Si). The abundance ratio of each element with respect to Fe were 
compared with those in stars with similar iron content of the so- 
lar neighbourhood. We confirm that O, Mg, Si, S and Na are con- 
siderably overabundant whereas Al, Ca, Ti and Ni appear to be 
practically consistent with the Galactic trends of these elements. 
We also report an upper limit of the Li abundance. 

These chemical abundances strongly suggest that the sec- 
ondary star captured part of the ejecta from a supernova ex- 
plosion that originated the black hole in Nova Sco 94. We have 
compared these element abundances with element yields from 
a variety of SN explosion models for different energies and 
geometries. An spherically symmetric explosion model of 30 
Mq progenitor (with a He core of 8.5 M©) suggests a mass- 
cut between 2-3.5 Mg for kinetic energies (1 - 20) x 10^' ergs, 
based on the relatively low abundances of Ca, Ti, Ni and Fe. 
The greatly enhanced abundances of Si and specially S in the 
secondary star favours higher kinetic energies and therefore a 
hypernova explosion, and requires efficient mixing processes be- 
tween the fallback matter and the ejecta. 

The kinematic properties of the system suggest that a na- 
tal kick was imparted to the compact object at birth due to an 
asymmetry in the neutrino emission if a neutron star formed 
first, and/or an asymmetric mass ejection. We have also in- 
spected a non-spherically synometric SN explosion model with a 
16 Mq He core progenitor, but this model provides unacceptable 
fits to the observed abundances because they require complete 
lateral mixing and produced too high abundances of Ca, Ti, Fe 
and Ni. We have also found relatively low Na and Al abundances 
which cannot be explained with the current spherical and aspher- 
ical SN models. 

Acknowledgements. We are grateful to Hideyuki Umeda, Ken'ichi Nomoto, and 
Nozomu Tominaga for kindly sending us their spherically symmetric explosion 
models and several programs for our model computations. We also thank Keiichi 
Maeda for providing us with his aspherical explosion models, and for helpful 
discussions. We would like to thank Martin A. Bitner and Edward L. Robinson 
for kindly sending us the code LinBrod. We are grateful to Tom Marsh for the 
use of the MOLLY analysis package. This work has made use of the VALD 
database and IRAF facilities. J. I. ackonowledges support from the EU con- 
tract MEXT-CT-2004-0 14265 (CIFIST). This work has been also funded by the 
Spanish Ministry project AYA2005-05149. 



Baumiiller, D., Butler, K., & Gehren, T. 1998, A&A, 338, 637 
Beer, M. E., & Podsiadlowski, Ph. 2002, MNRAS, 331, 351 
Bitner, M. A., & Robinson, E. L. 2006, AJ, 131, 1712 
Boesgaard, A. M., & Tripicco, M. J. 1986, ApJ, 302, L49 
Boesgaard, A. M., & Tripicco, M. J. 1987, ApJ, 313, 389 
Brandt, W. N., Podsiadlowski, Ph., & Sigurdsson, S. 1995, MNRAS, 277, L35 
Brown, G. E., Lee, C.-H., Wijers, R. A. M. J., Lee, H. K., Israelian, G., & Bethe, H. A, 2000, 
NewA, 5, 191 

Buxton, M., & Vennes, S. 2001, Publ. Astron. Soc. Aust., 18, 91 
Castelli, F., & Kurucz, R. L. 2003, lAU Symposium, 210, 20P 

Ecuvillon, A., Israelian, G., Santos, N. C., Mayor, M., Villar, V., & Bihain, G. 2004, A&A, 
426, 619 

Ecuvillon, A., Israelian, G., Santos, N. C., Shchukina, N. G., Mayor, M., & Rebolo, R. 2006, 
A&A, 445, 633 

Foellmi, C, Depagne, E., Dall, T. H., & Mirabel, I. R 2006, A&A, 457, 249 
Foellmi, C., Dall, T. H., & Depagne, E. 2007, A&A, 464, L61 

Galazutdinov, G. A., Musaev, E A., Krelowski, J., & Walker, G. A. H. 2000, PASP, 1 12, 648 
Gilli, G., Israelian, G., Ecuvillon, A., Santos, N. C, & Mayor, M. 2006, A&A, 449, 723. 
Gonzalez Hernandez, J. I., Rebolo, R., Israelian, G., Casares, J., Maeder, A., & Meynet, G. 
2004, Apl, 609, 988 

Gonz^ez Hernandez, 1. 1., Rebolo, R., Israelian, G., Casares, J., Maeda, K., Bonifacio, P., 

& Molaro, P 2005, ApJ, 630, 495 
Gonzlez Hemndez, J. L, Rebolo, R., Israelian, G., Harlaftis, E. T., Filippenko, A. V., & 

Chomock, R. 2006, ApJ, 644, L49 

Grevesse, N., Noels, A., & Sauval, A. J. 1996, The sixth annual October Astrophysics 

Conference, ASP Conf. Ser. 99, 1 17 
van der Hooft, E, Heemskerk, M. H. M., Alberts, E, & van Paradijs, J. 1998, A&A, 329, 

538 

Israelian, G., Rebolo, R., Basri, G., Casares, J., & Martin, E. L. 1999, Nature, 401, 142 

Kifonidis, K., Plewa, T., Janka, H.-Th., & MUer, E. 2000, A&A, 531, L123 

Kurucz, R. L. ATLAS9 Stellar Atmospheres Programs and 2 km s"' Grid. (CD-ROM, 

Smithsonian Astrophysical Observatory, Cambridge, 1993). 
Kurucz, R. L., Furenild, I., Bratxlt, J., & Testerman, L. 1984, Solar Flux Atlas from 296 to 

1300 nm, NOAO Atlas 1 (Cambridge: Harvard Univ. Press) 
Maeda, K., Nakamura, T., Nomoto, K., Mazzali, P. A., Patat, R, & Hachisu, I. 2002, ApJ, 

565, 405 

Marsh, T. R., Robinson, E. L., & Wood, J. H. 1994, MNRAS, 266, 137 
Mashonkina, L., Kom, A. J., & Przybilla, N. 2007, A&A, 461, 261 

Mirabel, L E, Mignani, R., Rodrigues, I., Combi, J. A., Rodriguez, L. R, & Guglielmetti, F, 

2002, A&A, 395, 595 
Orosz, J. A., & Bailyn, C. D. 1997, ApJ, 477, 876 
Pavlenko, Ya. V., & Magazzd, A. 1996, A&A, 311, 961 

Piskunov, N. E., Kupka, R, Ryabchikova, T. A., Weiss, W. W., & Jeffery, C. S. 1995, A&AS, 

112, 525 

Podsiadlowski, P, Nomoto, K., Maeda, K., Nakamura, T., Mazzali, P, & Schmidt, B. 2002, 
ApJ, 567, 491 

Sbordone, L., Bonifacio, P, Castelli, E, & Kurucz, R. L., 2004, MSAIS, 5, 93 
Shimanskaya, N. N., Mashonkina, L. L, & Sakhibullin, N. A. 2000, Astron. Rep., 44, 530 
Shchukina, N.G. 1987, Kinemat. Phys. Celest. Bodies, 3, 36 
Shchukina, N.G., & Trujillo Bueno, J. 2001, ApJ, 550, 970 

Shahbaz, T., van der Hooft, R, Casares, J., Charles, P. A., & van Paradijs, J. 1999, MNRAS, 

306, 89 

Shahbaz, T. 2003, MNRAS, 339, 1031 

Sneden, C. 1973, PhD Dissertation (Univ. of Texas at Austin) 
Takada-Hidai, M., Takeda, Y., Sato, S. et al. 2002, ApJ, 573, 614 
Takeda, Y, Hashimoto, O., Taguchi, H. et al. 2005, PASJ, 57, 751 
Takeda, Y, & Honda, S. 2005, PASJ, 57, 65 
Tominaga, N., Umeda, H. & Nomoto, K. 2007, ApJ, 660, 516 
Umeda, H.,& Nomoto, K. 2002, ApJ, 565, 385 
Umeda, H.,& Nomoto, K. 2005, ApJ, 619, 427 
Wedemeyer, S. 2001, A&A, 373, 998 

Zhang, S. N.,Wilson, C. A., Harmon, B. A., Rishman G. J., Wilson, R. B. et al. 1994, lAU 

Circ.,6046 

Zhao, G., Butler, K., & Gehren, T. 1998, A&A, 333, 219 
Zhao, G., & Gehren, T. 2000, A&A, 362, 1077 



References 

Abia, C, & Mashonkina, L. 2004, MNRAS, 350, 1127 

AUende Prieto, C, Barklem, P S., Lambert, D. L., & Cunha, K. 2004, A&A, 53, 181 
Al-Naimiy, H. M. 1978, Ap&SS, 420, 183 

Andrievsky, S. M., Chemyshova, L V., Usenko, I. A., Kovtyukh, V. V., Panchuk, V. E., & 

Galazutdinov, G. A. 1995, PASP, 107, 219 
Asplund, M. 2005, ARA&A. 43, 481 

Bailyn, C. D., Orosz, J. A., McClintock, J. E., & RemiUard, R. A. 1995, Nature, 378, 157 
Balachandran, S. 1991, Mem. Soc. Astron. Italiana, 62, 33 
Baumuller, D., & Gehren, T. 1996, A&A, 307, 961 
BaumuUer, D., & Gehren, T. 1997, A&A, 325, 1088 



